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ABSTRACT 


A Monte Carlo code is used to explore the full effects of the convection electric field 
on distributions of upflowing 0 + ions from the cusp/cleft ionosphere. Trajectories of in- 
dividual ions/neutrals are computed as they undergo multiple charge-exchange collisions. 
In the ion state, the trajectories are computed in realistic models of the magnetic field and 
the convection, corotation, and ambipolar electric fields. The effects of ion-ion collisions 
are included, and the trajectories are computed with and without simultaneous stochas- 
tic heating perpendicular to the magnetic field by a realistic model of broadband, low 
frequency waves. In the neutral state, ballistic trajectories in the gravitational field are 
computed. The initial conditions of the ions, in addition to ambipolar electric field and 
the number densities and temperatures of 0 + , H + , and electrons as a function of height in 
the cusp/cleft region were obtained from the results of Gombose and Killeen [1987], who 
used a hydrodynamic code to simulate the time-dependent frictional-heating effects in a 
magnetic tube during its motion though the convection throat. The distribution of the ion 
fluxes as a function of height are constructed from the case histories. 

The results show that during disturbed times, with a cross polar-cap potential equal to 
100 kV in the Heppner- Maynard convection-potential model, a high flux (> 10 8 cm -2 s -1 at 
1000 km) of low-energy 0 + ions flow outward from the ionosphere. At the altitude of 20,000 
km, these ions reach energies approaching 15 eV and their pitch-angle distributions become 
very narrow (half- width ~ 5°). If these ions enter the region in which the broadband, low- 
frequency wave field is present, they become heated; their energies increase and their 
pitch-angle distribution broadens during this process. Experimental observations of both 
these types of upflowing ions are discussed. 

A further result is that the very strong electric fields perpendicular to the magnetic 
field in cusp/cleft auroral structures that were reported by Sugiura et al. [1984], Sandholt 
[1988], and Sandholt et al. [1989a; 1989b] can directly accelerate 0 + ions from the iono- 
sphere, even in the absence of the heating and diffusion processes. An ion with an intial 
energy of 0.6 eV and pitch-angle of 45° at 2000 km is accelerated to 105 eV at 20,000 
km if it remains within magnetic field lines, assumed to be equipotentials, that produce a 
transverse electric field of 0.25 V/m at 300 km. 
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I. INTRODUCTION 


0 + ions, among other species, regularly flow outward from the ionosphere in region 
of the dayside cusp/cleft. The processes responsible for the outflow of these ions are of 
great interest because, as discussed by Cladis [1986], these ions enter the magnetosphere in 
the region where they may subsequently become appreciably accelerated by the convection 
electric field and transported to the near magnetotail. The pressure of these ions in the 
magnetotail, which was computed recently by Cladis and Francis [1992], suggests that 
these ions may trigger a substorm through the development of velocity shear instabilities. 
Moreover, during disturbed times, these ions often supply the major energy content of the 
ring current [Cladis and Francis, 1985; Horwitz, 1987; Delcourt et al. 1990]. 

Various wave-particle interactions [Chiu et al., 1989; Crew et al., 1990, and references 
therein] have been proposed to explain the ionic outflow from this region, and some of these 
interactions may indeed be responsible for the more energetic conics that are observed at 
high altitudes. However, the convection electric field, which is particularly strong in the 
convection throat, may be principally responsible for the outflow of these ions by acceler- 
ating them to higher altitudes where they may be further accelerated by other processes. 
Such a prominent role for the convection electric field is indicated by the correlation of 
the ionic outflux with indices of enhanced ionospheric convection [Lockwood et al, 1985; 
Yau and Lockwood, 1988; Yau et al. 1988]. Also, the escape flux is accompanied at 
lower altitudes (< 2 Re) with the so-called upwelling ions — ions of various species (He + , 
0 + , N + , NO + and 0 ++ , but mainly 0 + ) — that move upward with velocities of several 
km/s. This general upwelling of ionized atoms and molecules indicates the action of a 
strong ionospheric heating process. Moreover, hydrodynamic calculations of the frictional 
heating of ions in a magnetic tube convected through the throat have reproduced several 
features of the upwelling ions. The results of Gambosi and Killeen [1987] obtained with a 
time-dependent hydrodynamic code are particularly promising. 

In the research described here, an ion transport code is used to extend the hydro- 
dynamic calculations to high altitudes in the magnetosphere. Hydrodynamic codes do 
not properly handle non-Maxwellian plasmas or the acceleration and expansion of ion dis- 
tributions into the magnetosphere. Moreover, the transport code includes the following 
additional processes whereby the convection electric field affects the motion of the ions: 

(1) The outward acceleration of ions, along the magnetic field, by the Coriolis force evoked 
hy Ve , the E x B drift of the ions. Cladis [1986] has shown that at altitudes higher than 
a few REthis is the dominant force that acts on upflowing ions from the ionosphere. 

(2) Energy diffusion due to collisions of the ions with neutrals and ions in the presence 
of the convection electric field. The hydrodynamic codes compute this Joule (frictional) 
heating, which is very important at altitudes less than about 500 km. But the process 
is also important at higher altitudes where the anisotropy of the ions must be taken into 
account. 
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(3) Effects of multiple charge-exchange collisions. Owing to the divergence of the mag- 
netic field, ions can acquire high upward velocity components through multiple exchange 
collisions. The convection electric field, which on the average maintains the perpendicular 
velocity components of the ions at nearly the same level, effectively pumps energy into the 
parallel component. Intervening momentum-transfer collisions, together with the actions 
of the /jV||i? force and the ambipolar electric field enhance the outward velocity of the 
ions. 

(4) Wave-ion interactions. The convective flow in the throat region is susceptible to shear 
instabilities that result in magnetic and electric field turbulence [e.g., Kintner, 1976]. Ions 
can be accelerated stochastically by temporal and spatial fluctuations comparable with the 
ion gyroperiod and gyroradius, respectively, in the frame of the drifting ions. 

The ion-transport properties are constructed from case histories computed with a 
Monte-Carlo type code. The Monte Carlo code follows individual ions/neutrals as they 
undergo multiple (up to 4) charge-exchange collisions. In the ion state, the trajectories 
are computed in realistic models of the magnetic field and the convection, corotation, 
and ambipolar electric fields. Moreover, stochastic effects due to ion-ion collisions and 
wave-particle interactions are included. In the neutral state, ballistic trajectories in the 
gravitational field are computed. 

The results show that during disturbed times, with a cross polar-cap potential equal 
to 100 kV in the convection -potential model of Heppner and Maynard [1987], a high flux 
(> 10 8 cm -2 s _1 at 1000 km) of low-energy 0 + ions flow outward from the ionosphere. At 
the altitude of 20,000 km, these ions reach energies approaching 10 eV and their pitch- 
angle distributions become very narrow (half-width ~ 5°). Here, these ions are referred 
to as “unheated”, because they have not been heated by the intense broadband, low- 
frequency wave field that is often present on auroral field lines [Gurnett and Frank, 1977; 
Gurnett et al., 1984]. If these ions enter the region in which the waves are present, they 
become heated; their energies increase and their pitch-angle distribution broadens during 
this process. Both these types of upflowing ions have been observed. 

A further result is that the very strong electric fields perpendicular to the magnetic 
field in cusp/cleft auroral structures that were reported by Sugiura et al. [1984], Sandholt 
[1988], and Sandholt et al. [1989a; 1989b] can directly accelerate 0 + ions from the iono- 
sphere, even in the absence of the heating and diffusion processes. An ion with an intial 
energy of 0.6 eV and pitch-angle of 45° at 2000 km is accelerated to 105 eV at 20,000 
km if it remains within magnetic field lines, assumed to be equipotentials, that produce a 
transverse electric field of 0.25 V/m at 300 km. 
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II Transport Calculations 


1. Computer Codes 

The Monte Carlo code used in this research was principally constructed of codes we 
had previously developed to compute the transport of ions in the magnetosphere [Cladis 
and Francis, 1992, 1985; Cladis, 1986; Cladis, 1988] and in the ionosphere [Chiu et al., 
1988]. The former (magnetospheric) codes compute either the guiding-center motion or the 
full gyration motion of ions in prescribed models of the geomagnetic and geoelectric fields 
and contains options to compute simultaneously the stochastic acceleration of ions due 
to a specified distribution of transverse electric field variations. The latter (ionospheric) 
codes simulate ion heating and multiple charge-exchange effects in the auroral ionosphere. 
Trajectories of multiple charge-exchange events are computed by following atoms both in 
their charged state (guiding-center motion) and neutral state (ballistic motion) as they 
repeatedly gain and lose charge. A probabilistic test at each A t is conducted to deter- 
mine whether a charge exchange reaction had occurred. Periodic boundary conditions 
are assumed, i.e., particles that move across a boundary of the convecting magnetic tube 
are reinjected with the same velocity space characteristics at complementary configuration 
space locations in the tube. 

Additional codes were added to compute the effects of ion-neutral and ion-ion colli- 
sions. The ion-neutral collisions code accounts for elastic, inelastic, and charge-exchange 
collisions of 0 + on ambient 0 atoms; charge-exchange collisions of 0 + on ambient H atoms; 
elastic and inelastic collisions of energetic 0 (from 0 + on O charge-transfer events) on am- 
bient 0 atoms; and charge-exchange collisions of energetic O atoms with ambient 0 + and 
ambient H + . Here, ‘ambient’ refers to the initial distribution of a specie in the tube. Be- 
tween collisions the ion motion is assumed to be adiabatic; the ion is moved along and 
across the magnetic field by integrating the equations of the guiding- center motion de- 
scribed by Northrop [1963] (see Cladis [1988]). The collision probability is determined as 
described by Cladis and Francis [1985]. Since the ion-neutral collision time is much longer 
than the ion gyrofrequency, the velocity vector of an ion prior to a collision is determined 
by a uniform sampling (0° to 360°) of the ion gyrophase. The scattering angles appropri- 
ate for the particular collision are also randomly selected, and the calculation is resumed 
starting with the new velocity vector. 

Since the ion-ion collision frequency is very high, a stochastic process is used to de- 
termine the diffusion of the ions. This process is described in the following subsection. 

1.1 Ion-Ion Collisions 

The diffusion of 0 + due to ion-ion collisions is based on the diffusion coefficients de- 
rived by Chandrasekhar [1942] and Spitzer [1956]. Three coefficients are used: < (£va) >, 
< (£v s ) 2 >, and < (6v p ) 2 >. Here, < (Sva) >, the coefficient of dynamical friction, is the 
rate at which the test particles (the traced ions) are slowed down by collision with the field 
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particles (0 + ions of the ambient plasma); < (6v s ) z > is rate of increase of (Av) 2 in the di- 
rection parallel to the initial direction of motion of the test particles; and < (£v p ) 2 > is the 
rate of increase in the direction perpendicular to the initial direction of the test particles. 



< (Av,) >= _^G(V/V lh ) 

V th 

(1) 


< (Av,) ! >= ^G(V/V lh ) 

(2) 


< (Av p ) ! >= ^[*(V/V, h ) - G(V/V t h)] 

(3) 

where 

Ad = 8 ”* n(0 + )lnA 
nr 

(4) 

$(x) is 

the error function, 


$(x) = —i [ exp(— t 2 )dt 
V 7r Jo 

(5) 

and 

p/_x $(x)-x$'(x) 
{ > ~ 2x 2 

(6) 


where the prime denotes the derivat ive with respect to x. Additionally, V is the ion speed, 
V t h is the thermal speed (y / 2kT/m), e is the unit charge (4.8xlO -10 esu), and 


A = 


1.5, k 3 T 1 T 2 \* 

e 3 7rn e / 


(7) 


Here, Ti and T are the temperature of the field ions and the test-particle ions, respectively, 
and lie is the electron number density. In the region where the ion-ion collisions are 
important, T\ ~ T and n e ~ n(0 + ). 


As described above, the motion of an ion is assumed to be adiabatic during the time 
interval At between collisions. At the end of the time interval, the ion velocity vector 
in the Vg frame is altered by the velocity increments 6 \ d , £v s , and £v p , arising from the 
coefficients < (Av d ) >, < (Av s ) 2 >, and < (Av p ) 2 >, respectively. Here, 

6v d =< (Av d ) > At (8) 

and the velocity increments 6v s and £v p are determined probabilistically as follows: 

The distributions of Sv^, where the subscript k stands for s or p, are assumed to be 
Gaussian, 

p( * Vk) = W 1 (9) 
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where 


<j\ - 28i < (Av k ) 2 > (10) 

Then a random number Q, uniformly distributed between 0 and 1 , is selected and equated 
to the cummulative probability integral, 


Q = f P(y)dy 

J — oo 


where y = . After solving ( 11 ) for £, the velocity increment is given by 


( 11 ) 


<$v k = V2a^ 

Note that 6 v k may be positive or negative. 


( 12 ) 


The increments 6 vd and 6v s , which alter the magnitude of the velocity vector, alter 
the parallel and perpendicular velocity components, V|| and Vj_, by the amounts ( 6 v,j + 
6 v s )cosa and (£va + £v s )sina, respectively. The increment 6 v p changes both the pitch 
angle a and the guiding-center position perpendicular to B. Since the direction of 6v p is 
distributed uniformly in azimuth about V, the average rms value of <5v p along two orthog- 
onal axes normal to V is 8\ p /y/ 2 . Hence, by talcing one of these axes to be in the plane of 
V || and Vj_, it is clear that the pitch-angle change will be 


6a « 


1 6v p 

y/2 V 


(13) 


And the increment along the axis normal to the V||, Vj_ plane changes the guiding-center 
position by 


6 r = 


1 6v p 

y/2 n 


(14) 


where fi is the angular gyrofrequency of the ion. Since any direction of 6r about B is 
equally probable, it is taken to be at angle 7 with respect to the convection electric field 
E c , where 7 is randomly selected from the range 0 - 360°. 


The change in the circular energy (mVx 2 /2) of the ion due to the displacement of the 
guiding center along E c is 

£We = eE c 6 rcos 7 (15) 

Hence, the change in Vj_ due to 8v p is 


(*Vx)6v p = Vj_ 



6W E 

mV^/2 



(16) 
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Finally, by combining the effects of all the increments, the full alterations of Vy and 
V i| are given by the equations, 


6V i| = (£va + 6v s )cos(a + 6a) 


and 


tfVj. = (6v<j + 6v s )sin(a + 6a) + Vj. 

In the equations above Vj_ is limited to a small value > 0 


i + 45w-i 


mV^/2 


(17) 


(18) 


A restriction on the use of the probability distribution (9) is that <7 k be less than V . 
This condition can be assured if the time step At is taken to be sufficiently small such that 


At < 


V 2 

2 < (Av k ) 2 > 


(19) 


A similar procedure is used to simulate energy diffusion due to wave-particle inter- 
actions. The distribution of the perpendicular energy increments is taken to be Gaussian 
with cr = 2AtDj_j_, where Dj.j_ = e 2 wj_P„/(8m) is the diffusion coefficient and P„ is re- 
lated to the power spectral density of the waves. Here, too, there is a restriction on the 
magnitude of At, viz., 

At < 

[Cladis and Francis, 1985]. In the computer program, At is selected to satisfy the smallest 
of the limits given by (19) for k = s and k = p and (20). 


W_L 

2Dj_j_ 


( 20 ) 


2. Magnetospheric Models 

The magnetosphere is modeled by (1) the “short-tail” version of the magnetic field 
model of Tsyganenko [1987], as improved and coded by David Stern of NASA/Goddard, 
with the model parameters set for Kp > 5~ and tilt angle = 0°, and (2) the convection 
electric field model of Heppner and Maynard [1987] that is generally appropriate in the 
northern hemisphere for (B y )iMF > 0. The polar-cap potential drop, AtP, is taken to be 
100 kV, which is appropriate for Kp 5° and AE ~ 600 according to Table 1 of Heppner 
and Maynard [1987]. Therefore, the models used here are appropriate for more highly 
disturbed times. The vector sum of the corotation electric field and the convection electric 
field is mapped from the reference plane in the ionosphere, where the Heppner-Maynard 
convection potentials are specified, to the ion positions in the magnetosphere assuming 
the magnetic field lines to be equipotentials [Cladis and Francis, 1992]. For the neutral 
atmosphere, the model of Anderson and Francis [1966] is used. 

As discussed in the Introduction, Gombose and Killeen [1987] used a hydrodynamic 
code to simulate the time-dependent frictional-heating effects in a magnetic tube during 
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its motion though the convection throat. These authors also simulated the effects of 
precipitating electrons, which principally affected their electron distribution. The number 
densities and temperatures of 0 + , H + , and the electrons as a function of height in the 
cusp/cleft region were obtained from their results for times greater than 5 minutes. These 
properties are used to evaluate the diffusion coefficients (Equations (l)-(3)), the O~*"source 
used in the transport calculations, and the ambipolar electric field. The ambipolar electric 
field E a = — d^g/ds is evaluated by differentiating the equation, 


n e 


= n G 


e e* B /kT e 


( 21 ) 


assuming the electrons to have a Boltsmann distribution and the plasma to be quasineutral 
(n e as nj). 


The power-spectral density of the waves is taken to be P„ = (mV/m) 2 /Hz, 

where v is the wave frequency and j3 — 0.857. This is the electric-field wave power that 
was measured in the cusp/cleft region with the DE 1 satellite on March 7, 1984 [Andre et 
al., 1990]. It is assumed to be the same over the altitude range 2000 km to 20,000 km, as 
implied lay the measurements of Gumett et al. [1984]. Note that the wave power at the 
0 + gyrofrequency, which is used in the diffusion coefficient, increases with altitude as B - ^. 


3. Construction of Ion Properties 


The case histories of the traced ions are used to construct the following ion properties 
as a function of height: the directional flux j(w,a) in units of cm -2 s -1 sr -1 eV -1 , differential 
omnidirectional flux dJ(w)/dw in units of cm -2 s -1 eV -1 , integrated flux F up (cm -2 s -1 ) 
parallel to B, and number density n (cm -3 ). 


The ion flux j s (w,a) at the source altitude is taken to be Maxwellian with a tem- 
perature of 2400°. Each of the traced ions l mn represents a current of 6N mn ions per 
second moving upward along the magnetic field through a cross section £A S normal to the 
magnetic field. 


<5N 


mn 



27rsinacosada 


( 22 ) 


Here, the energy limits are chosen to divide the ion current into 5 equal parts between the 
thermal energy and infinity, and the pitch-angle limits divide the current into 5 equal parts 
between 90° and 180°. (The ions are traced in the northern polar cap so ions in this pitch- 
angle interval move upward). The energies w m and pitch-angles a n at the mid-current 
values within the above intervals are the initial values of the traced ions l mn . Twenty- 
five ions with the same initial conditions are traced. Since the number of different initial 
conditions is also 25, a total of 625 ions are traced. 


The ion case histories are tallied in planes normal to B to a maximum altitude of 20,000 
km. The listing for each ion that crosses a plane includes the trajectory number 1, time, 
geographic position, energy, pitch-angle, and the magnetic and electric field components. 
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In each reference plane, ions are tallied in five energy bins and five pitch-angle bins , 
which have widths that depend on energy- and pitch angle- ranges of the ions that cross 
the plane. The directional flux at the energy w k and pitch angle aj is then given by the 
summation 

\ _ Z^m,n gmnkj 

J Wk > a j (5A p 7r(cos 1 2 aj + i — cos 2 Qj_i )£w k 


(23) 


Here, g mn jk = 6Nmn when the energy and pitch angle of the ion l mn fall into the energy bin 
k and the pitch-angle bin j, respectively. Since the magnetic flux is conserved, the ratio 
of the areas 6A s /SA p is equated to the ratio of the average magnetic-field intensity in the 
plane to that at the source. 


The other properties of the ions follow directly from the flux: 

= X^( Wk ’ a j) 2T ( cosa j+i - cosa j-i) 

F up = ^j(w k ,aj)7r(cos 2 a j+ i - cos 2 * aj_i )£w k 

jtk 


and 


n up = ^7T (cos 7 a j+ ^ - cos 

j,k 


Vi) / 


dw 


__ 

w i *\J 2 "IX? j 771 




(24) 

(25) 


(26) 


Preliminary computer runs were made to determine an appropriate altitude for the 
source altitude. It was found that at altitudes less than about 750 km, the change exchange 
rate appreciably reduced the outflux of 0 + . Moreover, at altitudes up to about 1500 km, 
ion-ion collisions were very effective in keeping the ions thermalized. Therefore, the source 
altitude was taken to be 2000 km, and the ion properties were computed at the altitudes 
4,000, 8,000, 12,000, 16,000, and 20,000 km. 


Ill RESULTS 


1. 0 + Distribution in Absence of Wave-Particle Interactions 

First, the 0 + properties were computed for the case in which wave-particle interactions 

were not present. Table 1 shows the average directional flux of these ions in the energy 

and pitch-angle bins at each altitude. Note that the energy of the ions increases toward 
increasing altitudes, although the energy does not exceed 14.7 eV., and that the flux is 

high, above 10 8 ions cm -2 s _1 sr -1 eV _1 , even at the altitude of 20,000 km. In Figure 1 

the logarithm of the flux is depicted in Vj_, Vy coordinates. Here, the bins are in velocity 

intervals corresponding to the energy and pitch-angle intervals. The scale in the figure 

gives the correspondence of the logarithm of the flux to the grey-tone shading. This figure 



Tabic 1. Flux (cm J s *sr ^V” 1 ) of "unhealed 0 + ions in velocity, pitch-angle bins at 
altitudes 4,000 to 20,000 km. 


ALTITUDE - 4000 KM 


0.399- 2.594 
2.594- 3.647 
3.647- 4.458 
4.458— 5.142 
5.142- 3.746 


110.33-124.33 

6 . 74 61E + 05 
1 . 0119E+06 
O.OOOOE+OO 
O.OOOOE+OO 
0 . OOOOE+OO 


124.33-138 .11 

2 . 2265E+06 
1 . 94 02E + O6 
3.0615E+06 
1 .3916E + 06 
O.OOOOE+OO 


3 8 . 11-131 .90 

3. 5236E+06 
8. B089E+06 
7.3407E+06 
5 . 28 53E + 06 
5 . 8726E+05 


151 .90-165.68 

7 . 7 706E + 06 
1 .8404E+07 
1 . 9222E+07 
4 .0898E+06 
O.OOOOE+OO 


165 . 68-179 . 4 € 

4 . 1961E+07 
8. 0694E+07 
2 . 69 99 E + 07 
1 . 075 9E + 06 
0 . OOOOE+OO 


ALTITUDE - 8000 KM 


Pitch-angle* (deg) 


159.53-163.61 


3.931- 4.944 
4.944- 5.782 
5.782- 6.513 
6.513- 7.170 
7.170- 7.771 


O.OOOOE+OO 
0.0000E+00 
4 . 2682E+05 
8.5363E405 
8.5363E+05 


163.61-167.69 

O.OOOOE+OO 
3 . 1 99 IE +06 
9.0642E+06 
1 . 7062E+07 
3 * 1 991E406 


167 . 69-171.77 

1 . 4597E+06 
8.7585E+06 
4 . 1603E407 
1 . 4 597E407 
O.OOOOE+OO 


171.77-175.85 

0 . OOOOE+OO 
3 . 704 0E + 07 
6 . 4 522E + 07 
1 .0754E+07 
O.OOOOE+OO 


175 . 05-179 . 93 

1 . 3947E4 07 
2 . 4407E+08 
2 . 0223E+08 
6.9734E+06 
0 .OOOOEIOO 


ALTITUDE - 12000 KM 


_Pitch-an gie» (deg) 
"velociliei (km/»)~ 


4.233- 5.770 
5.770- 6.977 
6.977- 8.003 
8.003- 8.912 
8.912- 9.737 


169.22-171.37 

O.OOOOE+OO 
3. 9915E+03 
5 . 9073E+O6 
1.5966E+06 
7. 9830E+05 


171.37-173.53 

5.0925E+05 
2.54 62E + 06 
2 . 007 9E + 07 
6 , 1110E+06 
5.0925E+05 


173.53-175.69 

1 . 4175E406 
1 . 1340E+07 
4 . 6778E+07 
4 . 2525E+06 
7 . 0876E+05 


175.69-177.84 

1 . 17C2E406 
4 .5873E407 
3 . 6463E+07 
1 .5291E+07 
O.OOOOE+OO 


177 .84-180 .00 

7 . 02 6 6E 4 06 
2 . 0377E 4 08 
2 - 1783E + 0B 
1.7567E+07 
O.OOOOE+OO 


ALTITUDE - 16000 KM 


_Pitch-a ng)e» (deg) 
Velocitiei (kni/«) 


173. 22-174. 5B 


5.410- 7.188 
7.188- B.605 
8.605- 9.021 
9.821-10.901 
10.901-11.884 


3 . 6557E+05 
2 . 92 4 6E + 06 
2.5590E+06 
7.3115E+05 
3 . 6557E+05 


174.58-175.93 

1 .0742E+O6 
3 .7404E+O6 
8.4340E+06 
4 .2170E+06 
2.3420E+O6 


175.93-177.29 

1 .3087E+06 
5. B091E+O6 
1 . 8976E+07 
1 .7013E+07 
7.8521E+06 


177.29-178.64 

6.529BE+06 
1 .5236E407 
3 . 8090E+07 
4 .2443E+07 
1 . 3060E+07 


178 . 64-1B0 .00 

1 .3011E407 
4 . 8792E+07 
1.3662E+08 
2. 2770E408 
1.6264E+07 


JMlcWng lci (deg) 

Vdodtiei (km/«) 
7.010- 9.090 
9.090-10.771 
10.771-12.223 
12.223-13.519 
13.519-14.702 


173.33-176.26 

1 . 8170E+06 
2.4226E406 
6 . 0566E405 
6 . 0566E+05 
O.OOOOE+OO 


ALTITUDE - 

176.26-177.20 

3 . 4 992E + 06 
2 . 3328E + 0 6 
1.1664E+06 
1 .9440E+06 
2 . 721 6E+06 


20000 KH 
177.20-178.13 

5. 4374E+06 

4.34 99E + 06 
3 . 2624E+06 
6. 5249E+06 

1 . 35 93E + 07 


17B.13-n9.07 

1 .4490E+07 
1 . 1773E+07 
1 . 630 1 E 4 07 
3 . 1696E407 
3.094 1E+07 


179 .07-1B0 . 00 

4 . 88 B5E + 07 
5. 4316E407 
6 . 24 64 E 4 07 
9. 5053E407 
1 . 4394E+08 



clearly shows that the pitch-angle distribution of the ions become narrower toward higher 
altitudes as the ions become more field aligned. It also shows that at each altitude the 
highest flux is at the smallest pitch angle, but at the higher velocities the distribution 
tends toward a conic, especially at altitudes below 16,000 km. 

Histograms of the differential omnidirectional flux, dJ/dw, in the altitude planes are 
shown in Figure 2. Since this flux is obtained by integrating the directional flux over 
pitch-angle, its magnitude decreases with increasing energy. Note also that the peak of 
the energy spectrum moves toward the highest energy bin as the altitude increases. This 
effect is expected principally because the force Ve *de/dt, which accelerate ions along 
the magnetic field, has a component V||Ve * de/ds, which is proportional to the parallel 
velocity component of the ion. Here, e is a unit vector along the magnetic field and ds is 
an incremental distance along the local magnetic field line. The energy spectrum at 4000 
km is roughly characterized by a temperature of about 0.45 eV. 

Figure 3 shows the parallel flux (upper panel) and number density (lower panel) of 
the ions as a function of height. In statistical studies this flux is generally normalized to 
the reference altitude of 1000 km assuming the flux to be inversely proportional to B. At 
that altitude, the flux shown in Figure 3 would be 1.30 x 10 8 . The magnitude of this flux 
is in good agreement with the DE 1 EICS measurements of Onions in the energy interval 
0.01 * 1 keV during active times (3 < Kp < 5) in the sector 09 * 15 MLT [Yau et al., 1985], 

2. 0 + Distribution in Presence of Wave-Particle Interactions 

The particle trajectory calculations were then repeated with wave-particle interactions 
included. The wave field used for these calculations is described in Section 11,2. The 
average fluxes of these ions in the energy and pitch-angle bins at the altitudes 4000, 8000, 
12,000, 16,000, and 20,000 km are listed in Table 2. Comparison of these fluxes with 
those of Table 1 reveals, as expected, that the wave heating appreciably increases the 
energies and pitch angles of the ions. Histograms of the logarithm of the flux in shades 
of grey in Vy, Y ± coordinates are shown in Figure 4. These histograms clearly show that 
at the higher velocities the pitch-angle distributions are conics at all altitudes. Note also 
that, whereas the pitch-angle distribution of the unheated ions shown in Figure 1 rapidly 
becomes narrower toward higher altitudes, the pitch-angle distribution of the heated ions 
remains broad even at the highest altitude. 

Histograms of the differential omnidirectional flux, dJ/dw, at the same altitudes are 
shown in Figure 5. Here, note that the peak of the spectrum remaines in the lowest energy 
interval; it does not move toward higher energies as in the case of the unheated ions (see 
Figure 2). However, the spectrum becomes harder as the altitude increases. 

The parallel flux and the number density of the ions as a function of height are shown 
in Figure 6. Both these properties are similar to those of the unheated ions shown in Figure 

3. So these, too, are in agreement with the the DE 1 EICS measurements of 0 + ions in 
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Table 2. Flux (cm 2 6 1 sr*" 1 eV *) of “heated 0 + ions in velocity, pitch-angle bins at 
altitudes 4,000 to 20,000 km. 


'itch-wigles (deg) 


Velocities 

(km, 


3 

.07 - 

22. 

.52 

22 

.52 - 

31. 

61 

31 , 

.61 - 

38. 

61 

38 

.61 - 

44. 

53 

44 . 

.53 - 

49. 

75 


103. 96-110. 78 

2 . 5795E+05 
1 . 4 539E+05 
5 . 1589E+04 
9 . 3799E+03 
1. 4070E+04 


ALTITUDE - 

118.78-133.61 

2. 8717E+05 
1 .0018E+05 
5 . 0088E+04 
2 . 0035E+04 
1 .0018E+04 


4000 KM 

133,61-148.43 

1 . 9200E+05 
6.5O05E+O3 
0 . OOOOE+OO 
0. OOOOE+OO 
0 .OOOOE+OO 


148.43-163.26 

2 . 04 57E + 05 
0. 0000E+00 
0. 0000E+00 
0. OOOOE+OO 
0. 0000E+00 


163 . 26-178 . 08 

0 - 05 7 5E + 0 5 
0 .OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 





ALTITUDE - 

0000 KM 



Pitch- wig! 

les (deg) 









109.97-123.84 

123.04-137.71 

137.71-151.58 

151.58-165.45 

165.45-179.32 

Velocities 

(km/sT 






8.10 - 

35.49 

1 .6632E + 04 

2 . 6596E+04 

3 . 924 9E + 04 

7.2470E404 

2 . 73 95E * 05 

35.49 - 

49.53 

7.9831E+03 

4 . 8850E+03 

4 .5506E + 03 

5 . 51 40E + 03 

2.0444E+03 

49.53 - 

60.39 

5 . 3221E+03 

1 . 6283E+03 

1 .7065E+03 

7 . 0771E+O2 

0. OOOOE+OO 

60.39 - 

69.58 

1 . 9958E+03 

1 .0856E+03 

1 . 1377E+03 

0. OOOOE+OO 

0, OOOOE+OO 

69.58 - 

77.69 

0. OOOOE+OO 

1.0856E+03 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 


Pitch- wig] e* (deg) 
Velocities (km/ij 


10.40 - 50.47 
50.47 - 70.61 
70.61 - 86.16 
86.16 - 99.31 
99.31 - 110.91 


109.76-123.68 

2 . 5363E+03 
3 . 1704E+03 
1 . 1097E+03 
0. OOOOE+OO 
0. OOOOE+OO 


ALTITUDE - 

123.69-137.60 

5. 6682E+03 
9. 0176E+02 
7. 7293E+02 
1 .2882E402 
1 .2882E+02 


12000 KM 

137.60-151.52 

9 , 7038E+03 
9. 4343E+02 
4 , 0433E+02 
0. OOOOE+OO 
0. OOOOE+OO 


151.52-165.45 

1 . 5115E+04 
3 . 7320E+ 02 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 


165. 45-179.37 

7 . 4359E+04 
4 . 8600E+02 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 


Pitch-wigles (deg) 
Velocities {km/ij~ 


13.94 - 59.10 

59.10 - 82.41 

82.41 - 100.44 
100.44 - 115.70 
115.70 - 129.17 


108.28-122.50 

9.3798E+02 
4 .0199E+02 
2 . 6799E402 
1 . 3400E+02 
6 . 6998E+01 


ALTITUDE - 

122.50-136.71 

2. 4306E+03 
7.3975E+02 
4 . 2271E+02 
1.0569E+02 
5.2839B+01 


16000 KM 

136.71-150.93 

2 . 9413E+03 
2 . 1787E+02 
1 .6340E+02 
0. OOOOE+OO 
0. OOOOE+OO 


150.93-165.15 

6. 8091E+03 
2 . 9930E+02 
7 . 4B25E+01 
0. OOOOE+OO 
0. OOOOE+OO 


165.15-179.36 

3 . 03 1 9E+04 
7 . 7740E+02 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 


Pitch-wig les (deg) 
Velocities (km/s}~ 


19.21 - 62.17 

62.17 - 85.80 

85.80 - 104.20 
104.20 - 119.81 
119.81 - 133.60 


ALTITUDE - 20000 KM 


110.51-124.33 

4 . 5926E+02 
4 . 9753E+02 
1 . 91 36E+02 
2 . 2963E+02 
3 . 8271E+01 


124.33-138.15 

1 . 1362E+03 
5 . 6809E+02 
2 . 2092E+02 
3 . 1560E+01 
1 .8936E+02 


138.15-151.98 

1 . 7662E+03 
1 . 6663E+02 
9. 9976E+01 
6. 6651E+01 
0, OOOOE+OO 


151.98-165.80 

3 . 8179E+03 
9 . 31 10E+O1 
0. OOOOE+OO 
0. OOOOE+OO 
0 .OOOOE+OO 


165.80-179.62 

1 . 8147E+04 
1.1187E+03 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 



ALTITUDE =4000 KM ALTITUDE ^8000 KM ALTITUDE =12000 KM 
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Figure 1. Flux (cm 2 s 1 sr 1 eV *) of “unheated” 0 + ions shown in shades of grey in V 
V ii coordinates. 
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Figure 4. Flux (cm 2 s 1 sr 1 eV 1 ) of “heated” 0 + ions shown in shades of grey in V 
V n coordinates. 
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the energy interval 0.01 - 1 keV during active times (3 < Kp < 5) in the sector 09 - 15 
MLT [Yau et al., 1985]. 


IV SUMMARY AND DISCUSSION 

This study reveals that, during disturbed times, 0 + ions that are convected into 
the auroral zone may be accelerated along the magnetic field to velocities exceeding the 
gravitational escape velocity. The convection electric field, which is enhanced everywhere 
during disturbed times and especially in the convection throat, increases the outward 
flow through Joule heating in the altitude range 100 - 400 km and accelerates the ions 
along the magnetic field by the Coriolis force, Ve -de/dt, at all altitudes. The ions are 
further accelerated by the ambipolar electric field, which is enhanced principally by the 
precipitating electrons. The characteristics of the outflowing ions resulting from these 
effects are shown in Figures 1-3. Such upflowing-ion characteristics are similar those 
measured with the RIMS instrument on DE 1 [e.g., Lockwood et al., 1985; Moore et al. 
1986]. 

Within auroral flux tubes broadband, low frequency waves are generally present [Gur- 
nett and Frank, 1977, 1978], which accelerate ions perpendicular to the magnetic field. 
Thus, during the time interval that ions are convected through the auroral zone at appro- 
priate altitudes, they become heated by the waves. The power spectral density of these 
waves at the resonant frequency increases with altitude, so the ion-heating rate is higher 
if the ions are convected through the wave region at higher altitudes. Since the convection 
electric field model BC of Heppner and Maynard [1987] was used in the present study, the 
ions remained in the auroral region during the entire time (< 1500 s) required for them 
to be transported from 2000 km to 20,000 km. Hence, the ions were heated continuously 
during this period. The resulting ion distributions due to this wave heating, in addition 
to the effects described above, are shown in Figures 4-6. Here, the spectrum of the 
omnidirectional flux may be characterized by a “temperature” that increases from about 
21 eV at 4000 km to about 150 eV at 20,000 km. Moreover, the pitch-angle distribution 
remains broad at all altitudes to 20,000 km. It does not become narrower toward higher 
altitudes, as does the distribution of the unheated ions. Ion distributions displaying such 
“transversely-heated” distributions have been observed by many spacecraft at altitudes 
higher than about about 1 Re- Recently, a large data base of such distributions has been 
obtained with the Akebono spacecraft at altitudes less than about 5000 km. The SMS 
instrument on this spacecraft regularly observes transitions of upflowing ions from “un- 
heated” to “heated” distributions when the satellite enters the wave regions [Whalen et 
al., 1991]. 

Two sets of observations that substantiate these results are particularly noteworthy, 
those reported by Andre et al. [1990] and by Peterson et al. [1992]. In the former, both 
the “unheated” and “heated” 0 + ion populations were observed with the DE 1 satellite 
near 4 Re in the cusp/cleft region. The unheated ions had energies less than 65 eV and 
were highly aligned with the magnetic field in the region where the wave band was narrow 
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and the intensity was low. In the region where the wave field was broad-band and intense, 
the ion distributions were similar to the one shown in Figure 4 at 20,000 km. During these 
observations Ve was estimated to be 5 km/s at the DE 1 altitude. In our models the value 
of Vg at 20,000 km was 4.64 km/s. 

Peterson et al. [1992] describe simultaneous observations of upflowing ions at low and 
high altitudes near the same field line. One set of observations was made in the cusp/cleft 
region when Akebono was at an altitude of 5000 km and DE 1 was at an altitude of 
21,600 km. Ve was estimated to be 7 km/s at the DE 1 altitude. The SMS instrument on 
Akebono measured a maximum 0 + ion flux of 2 x 10 5 ions/(cm 2 s sr eV), with Tx ~ 10 
eV and T|| < 10 eV. The maximum directional flux at 4000 km, in the energy interval 
1.25 - 42 eV and pitch-angle range 163° - 178°, as listed in Table 2 is a factor of 4 higher 
than this value. Because of the viewing geometry of the SMS aperture, it is expected that 
the maximum flux was higher than the measured value. The perpendicular and parallel 
“temperatures” of the distributions are in satisfactory agreement. The EICS instrument 
on DE 1 measured maximum 0 + flux of 5 x 10 s ions/(cm 2 s sr eV) with Tx ~ 100 eV and 
T|| « 200 eV. The maximum directional flux at 20,000 km listed in Table 2, in the energy 
interval 30.7 - 322 eV and pitch-angle range 166° - 180°, is a factor of 3.6 higher than this 
value. The perpendicular “temperature” is approximately 100 eV, in agreement with the 
measurement, but the parallel “temperature” is about 140 eV, somewhat lower than the 
measured value. In view of the uncertainties in the model parameters used in this study, 
these measurements are considered to provide good support for the acceleration processes 
modeled in this study. 
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